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Abstract 

We consider motion in a periodic potential in a classical, quan- 
tum, and semiclassical context. Various results on the distribution of 
asymptotic velocities are proven. 

1 Introduction 

For a function V : M. d — > K which is periodic on a regular lattice C C M. d we 
study the evolutions 

ihd t W(t) = H h W(t) } W(0) = Id on L 2 (R d ), (1) 
d t ¥ = X H o $*, $° = Id on R 2d , (2) 

where H h = -f A+V(q), H(p, q) = ^ + V(q) and X H (p,q) = {-VV(q),p). 
This is done in the limits t — > oo and h — > 0. 

Why should one get interested in well-known things? In view of currently 
very active research on transport (anomalous or not) in condensed matter 
physics it is firstly desirable to put on a firm mathematical ground the folklore 
that 'motion in crystals is ballistic'; secondly one should try to obtain at least 
semiclassical information on quantities like the distribution of asymptotic 
velocities. 
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We investigate the asymptotic velocity lim^oo q(t)/t and the asymptotic 
behaviour of q 2 {t)/t a where a = 2 characterizes by definition ballistic, a = 1 
diffusive motion. 

Our results and the skeleton of the article are as follows: 

In Sect. 2 we show for a large class of V that the quantum motion is ballistic 
(Theorem |2.3| ). That class of potentials is not optimal in view of singularities. 
It includes, however, the Coulomb case. The modulus of the asymptotic 
velocity is bounded from above in a natural way (Corollary |2.4| ). 

In Sect. 3 we treat the classical motion in smooth potentials. In d — 2 
dimensions the motion is ballistic for high enough energies E; for d > 3 this 
is true for initial conditions outside a set of measure ~ 1/ \fE. There always 
exist fast orbits (of speed ~ VE), with a dense set of directions. This is even 
true for the ergodic case, where the asymptotic speed is zero with probability 
one, whereas almost all orbits are unbounded. (Theorem |3.1| ). 

The motion is never Anosov (Theorem |3.3| ). 

In particular this means that a gas of particles in a 'periodic' container 
is not Anosov if the interactions are smooth. So there may be small regions 
of regular motions and it seems unlikely that such a gas is ergodic. 

To the contrary the planar motion in periodic potentials with Coulombic 



(— 1/r type) singularities is known to be of Anosov type and diffusive, [|12 
In Sect. 4 we show that the distribution of asymptotic velocities -which are 
zero on a set of full measure- is dense in a disk of approximate radius \fE. 
(Theorem |4.3| ). 

Concerning semi classics we show in Sect. 5 that the quantum asymptotic 
velocities are always contained in a thickened convex hull of the classical ones 
for h small; they concentrate in measure inside the convex hull of the support 
of the classical probability distribution (Thm. |5.3| ). 

In particular in the classically ergodic case the positive speed of the quan- 
tal motion is only a quantum fluctuation vanishing in the semiclassical limit. 
The same is known also to be true for Coulombic potentials, see 



The above results are basically consequences of the Birkhoff Ergodic The- 
orem. This is interesting insofar our technique is likely to be applicable in dif- 
ferent semiclassical situations. On the other hand more specific information 
is needed to prove sharper results concerning the distribution of asymptotic 
velocities. 

This is done in Sect. 6, where we consider separable potentials. There 
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we have fast semiclassical convergence to the classical velocity distribution 
(Theorem |6.1| ). 

Acknowledgements. A.K. thanks the CPT in Marseille and the members 
of the IHES (Bures) for their hospitality during his visit in May 1997. 



2 Quantum Ballistic Motion 



Now we shall prove that the evolution of a quantum system in a rather 
general periodic medium is ballistic and that the asymptotic velocity exists. 
The latter is related to the band functions. 

It is known that for a certain class of singular potentials the spectrum of 



the Hamiltonian is absolutely continuous, see Thomas ||21|| , Reed-Simon JLS 



Knauf jnj. It is a folk conjecture that absolute continuity implies ballistic 
motion. Our proof in d dimensions is based on Bloch theory. 

We consider potentials V : W 1 —> M which are periodic w.r.t. a regular 
lattice CcR d 



V(q + i) = V{q) ( q eR d ,£eC), 

so that we may consider it function V : T — > 
M d /C, and calculate its Fourier transform TV : C* 
on the regularity of the potential are: 



on the unit cell T : = 
> C. Our assumptions 



(A q ): d = 2 
d = 3 
d>3 



V G L p (T) with p > 1, 

V G L 2 (T) and 

T{V) G l p (C*) with p<(d- I) /{d - 2). 



(A q ) implies that V is form small with respect to —A and even operator 
small for d ^ 2, see B. 



Denote D := —iftV, then 



D 2 
2 



is defined by its quadratic form with form domain Q(H h ) 
H x (R d )i for d ^ 2 the operator domain is D(H h ) = H 2 (R d ). 



Q(-A) 
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We denote by 

0(t) := W*(t)OW(t) 

(with the solution W(t) := exp(—iH h t / K) of (JT])) the Heisenberg time evolu- 
tion of an operator O. 

The symmetries of H h allow for a decomposition with respect to the group 
of lattice translations: Let C* be the dual lattice with unit cell T* and denote 
by 

U:L 2 (& d )^L*(T,^L*{T t dq)) = f L 2 (T, dq) dl 



\t*\ v ') y T » v y |T*| 

the unitary operator defined by extension from Schwarz space of 

UiJ>(k, q) = m) k {q) := £ e'^+^q + 1) (ip E S(R d )). 

eec 

The following facts are known in the literature and will be used below: 
Theorem 2.1 Let V satisfy (A q ). Then 

1. UH h lJ- x = f® H\k)^ with H h (k) = \{D + hk) 2 +V on L 2 (T), 
with form domain Q(H h (k)) = iJ 1 (T); 

2. k I— > H h (k) is a Type (B) analytic family; 

3. the spectrum of H h is absolutely continuous; 

4- H h {k) has compact resolvent, H h {k) = E^=i ^n(^)^i»W where E%(k) 
are the eigenvalues in ascending order, P^(k) the eigenprojections ; 

5. for every n the following are Lebesgue Nullsets: 
{k £ T* | E% is not differentiate at k}, 

{k E T* | P% is not differentiate at fc}, 

{k G T* | VkEftk) = 0}. 

Proof. (1-4) are proven in resp. in |H| for d = 2. (5) is proven in 
@, 0, see also §• □ 
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Remark 2.2 We emphasize that while the assumptions (A q ) are sufficient 
for absolute continuity, for d > 2 they are far from necessary for self- 
adjointness. It would be interesting to understand what happens in the gap! 

The result on ballistic transport in the quantum case is (see also the 
recent article || by Gerard and Nier): 

Theorem 2.3 Let V satisfy (A q ). It holds for ip with 
(Dip, Dip) + (qip, qip) < oo: 



Dip := lim 



q{t)ip 



t— »oo t 



/ J2p n (k)(D + hk)P n (k) 

^* n=l 

(«© oo 
J ^ n=\ 



dk 

]t*i 

dk 
T*| 



Uip 
Uip- 



and: 



t^oo t 2 



Proof. The map 



/ Yl \^kE n \ 2 \\p n um\\i H T } y£ > o. 

</T* n=1 I I 



V"-™/ D(t)iPdt-U- l \ Y^P n (k)(D + hk)P n (k) 

1 J V- 7t * n=l 



|T*| 



is uniformly bounded from if 1 — > L 2 , consequently it is sufficient to prove 
the assertion for ip in a dense set; indeed firstly by form smallness of V we 
have the estimate 

|| {D)ip\\ 2 = (ip, (1 + D 2 )ip) < Cl \(iP, HiP)\ + c 2 \\iP\\ 2 < c 3 (ip, (1 + D 2 )ip) 



so 



T 

secondly 



i A|| < ||<D)(fr + i)- 1/a ||||(ir + *) 1 ^|| <c\\ip\\ Hh 



u~ 



A* ^ IT* l 
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= / T\\P n (k)(D + hk)(H(k) +i)- l / 2 P n (k)U(H + i)V 2 iP(k) 

n=l 

< [ \\(D + hk)(H{k)+t)- 1/2 \\ 2 \\U(H + t) 1/2 i/j(k)\\ 2 ^ < const 
J T * 1 1*| 



2 dk 

L 2 (T) KpTj 



2 

1 ' 



Let ^ such that 



17V 



n=l 



|T*T 



Z7V- 



The set of these is dense in H 1 (M d ). Then 



r 



U[ - I D(t)$ dt 



N 



J2Pn(k)(D + hk)P n (k) 



71=1 



dk 
IT*! 



Z7V 



< 



oo,JV „y 

V - / exp (i(E m (k) - E n (k))t/h) dtP m {k){D + hk)P n {k)U^{k) 



2V 



- ^ P n (fc)(D + hk)P n (k)U^(k) 



n=l 



L 2 (T) 



IT*! 



n=l ^ J ° 



X 

m — 1 

(T -> oo) 



(E m (k)-E n (k))t/h 



dt{D + hk)P n (k)Uip(k) 



dk 
jf*j 



by dominated convergence, which is applicable because 
||P m (/c) X^n" • • II = O0-/T) for m n, almost all k, and is uniformly ma- 
jorized by 

const- sup \\P m (k)(D + hk)P n (k)Uip(k)\\ 2 L2(J) 

n=l,...,N 

which is summable with respect to m and k. 

In p| it was shown that for if) 6 H 1 ^) D D(\q\): 

q(T)ip = # + / D(t)i)dt. 
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It follows 



lim 



T 



lim — 



u- 1 



D(t)ipdt 



J2Pn(k)(D + hk)P n (k) 



n=l 



|T*| 



It remains to establish the identity 

P„(A;)(P> + hk)P n (k) = h- l W k E n (k)P n (k) 
for almost every fceT*. But this follows from: 

V fc P#P = E(V k P)P + EP(V k P) + P(V k H)P = EV k P + P(V k H)P, 
V k PHP = EV k P + PV k E, V k H = h(D + Hk), all valid in the quadratic 
form sense for almost every k by Properties [2.1| . 

This was the first assertion; the second one is a consequence thereof. The 
positivity is inferred from Theorem EOT □ 



As a corollary we get an estimate for the group velocity in one band: 
Corollary 2.4 For every n there is a set of full measure of k's such that 



h 



{\\1p\\=l,p n {k)1lj=lp} 



Proof. Let ip G L 2 (T), k such that E n ,P n are different iable at k. 
\h- 2 V k E n {k)\ 2 \\P n (k)n 2 = \\P n (k)(D + hk)P n (k)i;\\ 2 

< \\(D + hk)P n {k)n 2 

= 2{E n {k)\\P n {kW - (P n (k)<ip,VP n (k)iP)) 
implies the inequality. □ 



3 Classical Motion: Smooth Potentials 

The classical motion in a /^-periodic potential V on M. d is described by 
Hamilton's equations (|2|) on phase space P := T*M. d for H : P — > R, 
H(p,q) = |p 2 + V(q). If F G C 2 (IR d ,IR) (as we assume in this section), 
the flow : P — ► P exists uniquely for all times i G R. 
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We will analyze its restrictions <3>^ := $* f E to the energy shells 
Y, E := H-\E). 

Alternatively we study motion on the phase space P := T*T over the 
configuration torus (and mark corresponding objects with a hat). Using the 
phase space projection II : P — > P arising from the projection 7r : IR d — > 
T = IR d /£ of configuration spaces, we thus consider the flow : P — > P 
generated by the Hamiltonian function H : P — > R, H o II — H, and its 



is 



compact energy shells I1 E :— H 1 (E) with the restricted flows $^ := $ 

The Liouville measures A of the phase space regions if _1 ([V m i n , E\) E e 
are now finite, a fact which enables us to use notions of ergodic theory. 
The energy scales V mili : = inf qm dV(q), 



E 



Kncan := / V{q)dq/\T\ 
Jt 



and Vmax := sup (JgR dl / (g) of the dynamics will be used repeatedly. 

As a consequence of Birkhoff 's Ergodic Theorem for A-almost all x e P 

1 f T 

v ± (x ) := lim — / p(t,x )dt 

exist and are equal. In this case we set v := v ± , and otherwise v :— 0, thus 
defining the asymptotic velocity 

v : P -> M d 

which is a measurable phase space function. 

We denote its lift v o II : P — > M d to the original phase space P by the 
same symbol and thus have 

lim = f (XoJ 

A-almost everywhere. 

<E>* is called ballistic at x G P if ■y(x) 7^ (observe that by the above 
definition this implies existence and equality of v^ 1 ). 

We are particularly interested in the energy dependence of asymptotic 
velocity and thus introduce the energy-velocity map 

A := (H,v) : P -> R d+1 . (3) 
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A is measurable and generates an image measure v := A A 1 on 

Example. In the simplest case V = of free motion v is a smooth measure 
on the paraboloid 

A(P) = {{\v 2 ,v) | v G R d }. 

Like in the above example, in the general case v is invariant under (h, v) \— > 
(/i, — v ), since the motion is reversible. 

The equality H = \v 2 is in the general case replaced by the estimate 

\V(X)\ < yfi{H{x) - Knin). 

For regular values E of the energy one may consider the probability distri- 
bution of the asymptotic velocities v w.r.t. the normalized Liouville measure 
Xe on the energy shell By the above bound this is supported within a 
ball of radius \JE — V^m- 

Unlike in the above example in general it is not expected to depend weak- 
*-continuously on E, See Remark |3T2|.|l| below. 

Here are our results on classical ballistic motion for V G C 2 -potentials: 

Theorem 3.1 1. For d = 1 the motion is ballistic at x = (p,q) G P iff 
E := H(x) > V max , with asymptotic velocity 

sign(p) 

v(x) - 



l-^ l Q {2{E-V{q)))-\dq 

(I > being the period of C). 

2. For d > 1 and E > V max there exists a set Be C for which the 
motion is ballistic, whose directions 

{v{x)/\\v{x)\\ \ xeB E } 

are dense in S' d_1 , with moduli 

V2(E- V max ) ^ Mx)ll ^ y/2{E _ Vm . m)7 {xeBE y (4 ) 

V ''mean 



9 



3. For d = 2 and V G C 5 (M d , M) there exists a threshold E th > V max above 
which the flows $^ (E > E t h) are ballistic Xe -almost everywhere. 

Eth is given by the following condition. For E > Eth there are two 
geometrically different minimal tori Tf,T| C (by 'geometrically 
different' we mean: not related by time reversal symmetry X(p,q) : = 

(-P, q))- 

4- We assume here that V is 3d times continuously differentiable. Then 
for d > 2 there exist a threshold energy E t h > V max and for E > E t h 
subsets B E cT^e of measures 

\e(B e ) > 1 - VEth/E 
such that on Be the motion is ballistic. 

5. If the flow on the energy shell is ergodic w.r.t. Xe, then v — with 
\ E -probability one. However, if in addition E > V ma3L! the trajectories 
are unbounded with probability one: 



x e s E 



lim sup 

T 



pit, x )dt 



oo 



1. 



6. For d > 2 there are smooth C-periodic potentials V and energies E > 
Vmax whose energy shell contains a set of bounded orbits of positive 
measure: 



X, 



x e s E 



lim sup 

T 



p(t, x )dt 



< oo 



> 0. 



Example. Consider first in d = 1 dimensions the potential V(q) = cos(g). 



With the formula of Thm. |37i] for E > 1 the asymptotic speed equals 
v(E) 



7T 



Ve~^T 



y/2EllipticJC(2/(l- E)Y 

(EllipticJC being the complete elliptic integral of the first kind) and v(E) = 
for -1<E<1. 

For d = 2 this leads to a distribution of asymptotic velocities for energy 
E of the potential V(q) = cos(gi) + cos(g2) depicted in Figure [I]. Observe 
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Figure 1: Distribution of asymptotic velocities v for the separable potential 
V(q) = cos(gi) + cos(g 2 ) and energy E = 3 

that there is a positive probability for motion along the axes. 



Proof. 

1. v(xq) = l/T where T := Jjj ^ i s the time needed for the spatial period 

2. The idea is to construct periodic but non-contractible orbits on the 
torus. These are covered by ballistic orbits in the configuration space 

K 

For E > Vmax we consider the geodesic motion on T in the Jacobi metric 



9e( 



:= (E-V(q))-J2dq t 



(5) 



i=l 



The geodesies of that metric are known to coincide with the solution 
curves t h- > q(t,Xo), x G Eg, up to a time reparametrization r i— > t(r). 

In every nontrivial homotopy class / of the fundamental group 7Ti(T) = 
£ we find a shortest closed geodesic c : S 1 — > T (with S* 1 := R/Z). The 
length 



L{c) = 



f 


dc(r) 


/o 


dr 



■ ^E - V{c{r))dT 



(6) 
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(||-|| denoting Euclidean norm) of that geodesic in the Jacobi metric is 
the infimum of the lengths of the curves in its homotopy class I. Let 
the corresponding solution curve 1 1— > q(t,x ) on the torus have period 
T. 

We have q(nT, xq) = q(0, Xo)+n-l for motion in configuration space R^, 
starting from a point xq G n _1 (xo). Therefore the asymptotic velocity 
v(xq) of this orbit exists and equals v(xo) = l/T. So our task is to 
estimate T from below and above. 

The upper bound 

llt?(x)ii < vw 



Vrr 



follows from the general bound < >/2(E — V min ) if (p, q) G S B . 

In order to prove the lower bound for ||u||, we derive an upper bound 
for the period T and argue as follows. The length L(c) of our minimal 
geodesic c is shorter than the lengths of all the homotopic straight lines 
c = dq : S 1 — > T, 

c(r) := q + t ■ I (mod C) 

starting from a point g 6 T. But by formula (|6]) the length L(c) = 
\\l\\- j'l - V{c{r))dr. By concavity of x i-> v 7 ^" 



/ y/E - V{c{r))dr < ^ E - V m 
Jo 

for some q. So we obtain 

L(C) < L(5) < |/| ■ V 7 ^ " K.ean- 

On the other hand the period 



(7) 



T 



lo dT dT 


r 






Jo 





/>/2(£7-K(c(r)))dr 



< L(c)/ (V2(E — Vmax)). 
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Together with (0) this gives T < \\l\\ ■ y/E — V meim / (V^(E — V max ) ) from 
which the lower estimate for the asymptotic speed ||u(xo)|| = ||Z||/T in 
(D follows. 

The asymptotic direction v(xo)/\\v(xo)\\ of our ballistic orbit equals 
l/\\l\\- But the directions of lattice points I 6 C \ {0}, seen from the 
origin are dense in S 1 . 

3. We show that under the existence assumption for the tori Ti, T 2 C 
the motion is ballistic A#-a.e. 

For E > V^ax and d = 2 the energy shell is diffeomorphic to S 1 x T 2 , 
S 1 representing the circle of directions p/||p||. 

By the minimality condition the tori 

Ti,T 2 C t E C P ~ M 2 x T 

(and their time inverses) X(Ti), X(T 2 )) project diffeomorphically to 
the configuration torus T. Thus we may represent them as graphs of 
functions Pi, P<i : T — > R 2 . For i — 1,2 the mapping 

A/IIPII :T^S X 

are the local direction, and is by minimality topologically trivial. So 
their complement 

S B -(TiUT 2 UJ(Ti)UX(T 2 )) 

in the energy shell of four components diffeomorphic to thickened two- 
tori. These components roughly correspond to sectors of directions in 
which the particle is forced to move. 

The problem which has to be overcome is that these sectors of directions 
depend on the point q, and that their union for all q may have a total 
opening angle of more than tt. Thus it may happen that the particle 
goes backward for some time. 

Without loss of generality we consider the component C C Tie — (T 2 U 
T 2 ) which consists of points (p, g) e which are linear combinations 
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p = a\Pi{q) + 0:2^2(9) of the points (Pj(g),g) € TTf with positive coef- 
ficients Oti. 

On the invariant tori T 2 the motion is conditionally periodic with fre- 
quency vectors u>i £ M 2 \ {0}. If we consider the Lagrangian manifolds 
Mi C Eg which under II project to the tori T 2 , these manifolds are 
not only diffeomorphic to R 2 , but by our minimality assumption they 
project under II : P — > IR 2 , (p, q) 1— > g diffeomorphically to the config- 
uration plane M 2 . Thus they induce two flows 

^ : r2^ r2) ^^nofojny- 1 . 

These are nearly linear in the sense 

y t i (q ) = q + u l -t + O(t ) (8) 

with 0(t°) uniform in go, since they come from a conditionally periodic 
motion on a torus. The flow lines of and ^2 both foliate the config- 
uration plane and are transversal to each other. Since both foliations 
project under n : IR 2 — > T 2 to foliations of the (compact!) torus, the 
angles under which these foliations intersect are bounded away from 
and 7r, see Figure |2|. However, and ^2 do not commute in general. 
Nevertheless, we may use them to find an adapted coordinate system 

* : M 2 -> M 2 , s 2 ) h-> o ^ 2 (g ) n ttf o ^(g ), 

since the orbit tyf o ^^(Qo) through the point \&2 2 (<7o) has a unique 
intersection with the orbit ^Z^ 1 o tyf (g ) through (go). By the above 
remarks we have 

s 2 ) = go + ^1 ■ si + W2 ■ s 2 + 0(1), 

and the Jacobian of \I/ is uniformly bounded. 

We consider the component C = n~ 1 (C') of and an initial point 
= (PO)9o) ^ C- By compactness of T the angle between P\(q) and 
P 2 (g) is bounded away from and 71". Thus for some c > each point 
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<h 

Figure 2: The two foliations of the configuration torus 

(p, q) € C has a momentum vector p which is a linear combination 
p = a\px + OL2V2 with (pi, q) e Mi and a\ + «2 > c. 

Thus the ^/-coordinates Si, S2 are increasing along the trajectory t 1— > 
q(t, xo), and there exists a c' > with 4:{si(t) + S2(t)) > c'. 

The linear coordinate g : M 2 — >■ R, g := (cl>i/||u;i|| + k^/H^H) ■ q on the 
configuration plane increases (at least) linearly along the trajectory. 
Namely, by (§) the trajectory meets the inequality 

9(s(Mo)) 

V 11^1 II 11^2 

= (|M| + ^) »i(t) + (|M| + ^) s 2 (t) + (9(1) 
> c//(si(i) + s 2 (i)) > c/c// • t, 

if i is large. That is, if the asymptotic velocities ^(xo) exist and are 
equal, they must be non-zero. 



) ■y(s l (t),s 2 (t) 
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The existence of such minimal tori for large E follows under our dif- 
ferentiability assumption V G C 5 (1R 2 ,M) from the results of the paper 
03 of Poschel. 



4. In the Jacobi metric (||) the perturbation of the integrable part equals 
(1 - V(q)/E) ■ Y?i=i d( li ® Mi- Since V G C 3d (T,M), the norm of the 
perturbation is finite and proportional to E^ 1 . The unperturbed part 
of the Hamiltonian function for geodesic motion in that metric is just 
the Hamiltonian of free motion. So up to a linear transformation the 
momenta coincide with the action variables of this integrable system, 
and the nondegeneracy condition of the frequencies is satisfied. These 
depend analytically on the action variables, and the perturbation is 
C 3d and of size e = 0(1/ E) We can apply the result by Poschel on the 



measure of KAM tori for a perturbation of order e (Corr. 2 of |17|), 
which says that the measure of the complement of the KAM tori is of 
order ■s/e. 

5. If \e({v = 0}) < 1, there exists an index j G {1, . . . , d} with probabil- 
ities 

X E ({vj > 0}) = \ E ({v 3 < 0}) > 0. 

But this contradicts ergodicity, since these two exclusive events are 
flow-invariant. 

To show that the trajectories are unbounded with probability one if 
E > Vmax, we consider the flow-invariant measurable events E n C 
n£N defined by 



E n := <x Q G E e 



3ti,t 2 G 



*2 

p(t, x )dt 



> n 



We know from part 2 of the theorem that there are ballistic trajectories. 
Each point x G on such a ballistic trajectory is contained in all the 
sets E n . By absolute continuity of || f f 2 p(t, x )dt\\ w.r.t. the initial 
condition xq we conclude that for all n G N the Liouville measure 
\E(E n ) > 0. But the flow being ergodic, and E n flow-invariant, A^(i? ri ) 
can only zero or one. So the set n ne N-En of unbounded trajectories has 
measure one. 



16 



6. To construct potentials V G C°°(IR d , M), d > 2, which have energies 
E > Knax with many bounded orbits on the energy shell, one writes 
V (<l) : = J2ieC W (l- £ ) with W G Cg°(R d ,R), W(q) := W(\q\), with 
W(r) = for r > |min£ e £ ||£||, so that the supports of the lattice- 
translated W do not overlap. So as long as the particle is captured near 
a lattice point £, the motion is one in a potential centrally symmetric 
around £, and the angular momentum around that point is constant. 

Thus one reduces the dimension by considering the effective potential 
Wi,(r) := W(r) + for angular momentum L. For a given choice 
of L / one chooses W < so that the effective potential has a 
strictly positive nondegenerate minimum r : Wl(tq) > 0, 4zWl{tq) = 
0, ^Wi(ro) > 0. Then the assertion holds for E = Wl(to) + e, since 

Knin < 0. 

Whether one can construct for d > 3 potentials with a positive measure 
of bounded orbits on energy shells of arbitrarily large energy, is a much 
more complicated question. □ 

Part 2 and 3 of the above theorem show that the minimal KAM tori play 
an important role in the distribution of asymptotic velocities. However, it is 
known [II] that for d > 2 an energy shell can only be foliated by such tori if 
V is constant. This suggests that tori which do not diffeomorphically project 
to the configuration torus are important for the high energy distribution of 
asymptotic velocity, see Figures and |[ Indeed, in the separable case 
V(q) = with Vj non-constant, the probability to move in each of 

the directions £i, . . . ,£d is positive (Figure |l|). In the non-separable case like 
the one depicted in Figure it can be proven by perturbation arguments 
that there are lattice-rational directions in which the particle moves with 
positive probability. If for every probability measure ue which is absolutely 
continuous w.r.t. A# 

D{v E ) := hm 

exist and is positive, we call $^ diffusive. 

Moreover, we say that the flow $^ is diffusive in the strong sense if 
q(t, xq) I y\t\ converges weakly to a Gaussian distribution with positive co- 
variance matrix. An example of a strongly diffusive flow is given in |L2] . 
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Figure 3: Motion on unperturbed (left) and perturbed (right) invariant torus 
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Figure 4: Motion on an invariant torus which does not project diffeomorphi- 
cally onto the configuration torus 
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Figure 5: Distribution of asymptotic velocities v for the non-separable po- 
tential V(q) = cos(gi) + cos(gi + q 2 ) and energy E = 3 (numerical) 

Remarks 3.2 1. Conversely to the third statement, ford = 2 a lowering 
of the energy may lead to the destruction of the second to last KAM 
torus, which in turn may lead to a discontinuous decrease of the group 
velocities. 

As only for d = 2 the d-dimensional KAM tori have codimension one 
in the 2d — 1-dimensional energy shell for d > 3 Arnold diffusion 
may lead to initial conditions of positive measure which are not ballistic. 

2. By the statement 4) there always exists a threshold energy above which 
the motion is not diffusive. 

3. The mean classical velocity v for d = 1 equals the velocity expectation 

{ip - iVip) 



Theorem 3.3 If d > 2, then there is no energy E for which $g is an Anosov 
flow. 



of the WKB function 
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Proof. If V = 0, the motion is integrable and thus never Anosov. So we 
may assume V non-constant. The Hamiltonian is optical, that is, strictly 
convex on each fibre. 

Then for E < V max the energy shell touches the zero section of T*M. d . 
Thus by Theorem 1 of Paternain and Paternain JT^| $^ is not Anosov. 



For E > Vmax Theorem 3 of |14| which generalizes a Theorem of E. Hopf 
says that the flow <3>^ has conjugate points if V is non-constant. Thus by 



Theorem 1 of 16 the flow cannot be Anosov either. □ 



Remarks 3.4 1. Of course Theorem \3. 5| does not imply that motions in 
smooth potentials on T cannot be ergodic. To the contrary, Donnay 
and Liverani gave in ^ a method to construct such ergodic C°°(T) 
potentials for d = 2 freedoms. These potentials were constructed in 
such a way that for a given energy they contained circularly symmetric 
pits with a parabolic circular trajectory. Any decrease of the energy 
then makes this trajectory elliptic and the motion non- ergodic. 

Our theorem shows that non-hyperbolic trajectories like in that example 
must necessarily appear. In general we conjecture that for d = 2 these 
trajectories lead to anomalous diffusion effects and are incompatible 
with diffusivity in the strong sense. 

Sinai and Kubo gave examples of repelling continuous potentials on a 
torus which lead to ergodic flows. However, in this case the potentials 
could not be chosen to be C 1 so that they, too cannot serve as coun- 
terexamples to our theorem. See /@/ for a discussion. 

2. Motion of k particles on an d-dimensional configuration space with 
periodic boundary conditions and mutual forces of potential type can be 
described by the motion of one particle on a k ■ d-dimensional torus. 
Thus Theorem 2 implies that it will be very hard to show ergodicity of 
gases if the interparticle forces are smooth. 

3. A geometric version of the above theorem is: Geodesic flows on a torus 
(T, g) are never Anosov. This follows from the generalizations of Hopf 's 
Theorem by Burago and Ivanov ||/ ; together with Theorem 1 of fJ^/ . 

This is clearly not a mere consequence of the topology of the unit tangent 
bundle S d ~ l x T, since for d = 2 this is a three-torus, and the simplest 
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example of an Anosov flow (a suspension of Arnold's cat map) is one 
on T 3 . 

4 Classical Motion: Coulombic Potentials in 

d = 2 

We now treat motion in a planar crystal with attracting Coulombic forces. 
We fix the locations of the nuclei within the crystal by selecting m > 1 points 
Si, . . . , s m e T> in the fundamental domain 

D : = { Xx t x + x 2 £ 2 | xi, x 2 e [0, 1)} C 

of the lattice £ctf with basis £±,£2- The nuclei attract the electron with 
the charges Z±, . . . , Z m > 0. That is, we assume the potential of the form 
V(q) ~ —Zi/\q — Si\ for q near Sj. Now by the periodicity of the crystal the 
potential is singular at the points of 

S :={si + £\ie{l,...,m},£eC} 

and thus only defined in the punctured configuration plane M := M.^\S. 
Sometimes we identify the plane with C. 

Definition 4.1 A potential V G C°°(M, R) which is C-periodic 

V(q + £)=V(q), ( q eM,£eC) 
is called Coulombic if for e > small the functions fi, . . . , f m , 

f(n \ - I V( Si + Q 2 )-QQ ,0<|Q|<5 

are C°°. 

The reason for the somewhat odd-looking definition is that we want to regu- 
larize the Coulomb singularities by using the so-called Levi-Civita transfor- 
mation. 

Observe that V max := sup 96A ^ V(q) and Kncan ar e still well-defined finite 
quantities. 
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The classical motion is generated by the Hamiltonian function 

H(p,q):=±p 2 + V(q), {(p,q)eT*M). 

Due to collisions with the singularities in S the Hamiltonian flow on the 
cotangent bundle T*M of the punctured plane M does not exist for all times. 

However, as described in Lemma [4.2| below, the flow can be smoothly 
regularized. 

Lemma 4.2 There exists a unique smooth extension (P,u,H) of the Hamil- 
tonian system (T*M, dqi A dpi + dq 2 A dp%, H), where the phase space P is a 
smooth four-dimensional manifold with 

P := T*M U [J R x S 1 (9) 
s 

as a set, u is a smooth symplectic two-form on P with 
uj \ t *m = dqi A dpi + dq 2 A dp 2 , 

and H : P —>■ R is a smooth Hamiltonian function with H \ t *m = H . 
The smooth Hamiltonian flow 

: P -»• P, (t E R) (10) 

generated by H is complete. 

For all energies E which are regular values of V, the energy shell 

E E := {x G P | H(x) = E} (11) 

is a smooth, three-dimensional manifold, and we write $^ := $* ^ s . 

Proof. The construction works locally near the singularities s G S. We 
shortly explain the method by considering the simplest case. For more details 
see Prop. 2.3 of 0, where a scattering potential is considered. 
One may linearize the Kepler flow with Hamiltonian function 

H K {p,q) :=\p 2 -Z/\q\ 

in a suitable neighbourhood of R^ x {0} in the phase space T*(R^ \ {0}), 
using the canonical coordinates Hk, T, L and <p, where L(p,q) := p A q is 
the angular momentum, T(p, q) is the time needed to come from the phase 



22 



space point (p, q) to the pericenter of the Keplerian conic section and <p(p, q) 
is the angle between the direction of that pericenter and, say, the 1-direction. 
Except T these phase space functions are constant along the Kepler flow, and 
the collision points correspond to T = 0, L = 0. The remaining coordinates 
(Hk, (p) take values in a cylinder K x S . 

Similarly, because the singularities of V are of the Kepler form, one may 
thus complete the phase space T*M by gluing one cylinder for each singu- 
larity in S. □ 

So after this regularization the Hamiltonian flow exists for all time, and we 
are in a similar situation as in the case of a smooth periodic potential treated 
in Sect. |3|. In particular we also consider the motion over the configuration 
torus T. 

Theorem 4.3 If A\n(E th — V) > 0, then for all E > E th the intersection 
of the set C M 2 of asymptotic velocities for energy E with the disk of 

radius 

V2(E - Knax) 
yE V^nean 

is dense. 



Remark 4.4 Since the motion is diffusive ^T^J and thus in particular the 
asymptotic velocity can only be non-zero on a set of measure zero, this should 
be considered as a 'very large deviation' result. 

The ballistic orbits which we construct are of 'stop and go' type, that is, 
they periodically change between fast motion in a given direction and localized 
motion. Such orbits are somewhat special to our Coulombic potential and do 
not exist in the general case. Compare, however, with Thm. \3. 4 H- 



Proof. Although Lemma [4.2| solves the problem of collision orbits and reg- 
ularizes the motion without changing it otherwise, we will now introduce a 
second regularization which is more useful when one tries to construct orbits 
with prescribed asymptotic velocity, whereas it leads to a different phase 
space and a new time parametrization. 

For Coulombic potentials, the Jacobi metric (||) becomes singular near 
the positions S of the nuclei. However this is only a coordinate singularity, 
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and Gaussian curvature does not diverge. This can be seen by considering 
the case of the Kepler flow and the regularizing complex coordinate Q G C 
with Q 2 = q, that is, the Levi-Civita transformation. It turns out that in 
this coordinate the Jacobi metric, originally defined only for Q ^ 0, can 
be smoothly extended to a non-singular Riemannian metric on the whole 
complex Q plane. Note that the Q plane is a two-fold branched covering 
of the original q plane, with a branch point at the position q = of the 
singularity. 

In |12| this local construction was globalized using the toral Riemann sur- 
faces T = M/C, its four-fold covering torus T 2 := M/(2C) (with projection 
11x2 j : T 2 — > T) and the compact Riemann surface 

n 2 _ UT=i <Kg - gj + iiMg -Si + h + 1 2 ) \ n9 x 



M 4 := <^ {q, Q) G T 2 x P 



:= 2 ] [ (1 — z/w) exp (z/iu + ^(z/w) 2 ) 

w&2C\{0} 

being the Weierstrass cr-function and P := CUoo the Riemann sphere. The 
map II M4jT 2 : M 4 — > T 2 , (q,Q) i— >■ g is a two-fold branched covering with 
branch points at the singularities, all branch numbers equalling one. 

Thus by the Riemann-Hurwitz relation the genus G(M±) of M4 equals 
Q{M^) = 2m + 1 (m being the number of singularities in the fundamental 
domain). Since thus the genus is > 3, by Gauss-Bonnet the integrated Gaus- 
sian curvature § M ^K E dM/± = —4rr ■ (£?(M 4 ) — 1) of the lifted Jacobi metric 
g^E on M 4 becomes negative. Due to the branched covering construction 
the metric g^E, originally not defined at the branch points, can be smoothly 
extended to these points by taking limits. So the geodesic flow 4>\ E on the 
unit tangent bundle TxM 4 of the surface (M 4 , g^) is smooth and defined for 
all times. 

In terms of the potential V the Gaussian curvature Ke equals 

yn _ g: V(g))AV(g) + (WQ?)) 2 A\n(E th -V(q)) 

KE[q) ~ 2{E-V{q)f " 2(E-V(q)) (13) 

with A and V denoting the Euclidean Laplacian and gradient, respectively. 

For many Coulombic potentials V the Gaussian curvature of M 4 becomes 
strictly negative (K E (x) < for all x G M 4 ) if E is large enough. Clearly this 
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can only happen if AV < in (|I3"D, and negativity of ( p~3| ) is then preserved 
if one enlarges the energy E. 

Similar to the proof of Thm. 3TT|j2] we construct ballistic orbits in IR 2 by 



finding closed geodesies of minimal length. However, due to the presence of 
singularities, these geodesies are not constructed on T but on the smooth 
Riemannian surface (M 4 , g^)- 

Since (M 4 , g^E) is closed, there exists a closed geodesic in every non-trivial 
conjugacy class of the fundamental group 7Ti(M 4 ). Moreover, if the Gaussian 
curvature Ke is strictly negative (as it is the case if A ln^th — V) > and 
E > Eth), then this geodesic c : S 1 — > M 4 is essentially unique within its 
conjugacy class. Namely, any other closed geodesic in that class coincides 
with c up to a shift of the initial point (see, e.g., Klingenberg (TT|, Thm. 
3.8.14). 

As we are interested in ballistic trajectories on M 2 , we seek closed trajec- 
tories on T which are, however, non-contractible curves on the configuration 
torus. How many such orbits can we construct by projecting closed geodesies 
in M 4 to T? 

This question can be answered by considering the covering projection 
II .!/;.• := n T 2 jT o n M4jT 2. (14) 
This continuous map induces a homomorphism 
(n M4 , T ) t : 7Ti(M 4 ) -> 7Ti(T) 

of fundamental groups. We claim that the image subgroup Im C ~K\ (T) equals 

Im = 7T!(T 2 ) = 2£, (15) 

that is, consists of the equivalence classes of all loops c : S 1 —>■ T which 
surround the torus T = S 1 x S 1 in both basic directions an even number of 
times (this statement is of course independent of the base). 

([15]) follows from the definition ([14]) if we show that the homomorphism 

(n^), : tt x (M 4 ) -> 7n(T 2 ) 

is onto. To show this, we consider loops in T 2 based at b G T 2 , where we 
assume the base point b not to be a one of the 4m branch points of n M4jT 2, 
see (|T2"P, so that it has exactly two preimage points 60, &i £ Mi. 
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Now we can uniquely lift any based loop c : S 1 — > T 2 which avoids the 
branch points, to a path c : [0, 1] — > M 4 starting at c(0) := bo and ending 
either at 5(1) = b or at c(l) = &]_. 

On the other hand, we can find one fixed based loop I : S 1 — > T 2 which 
is contractible in T 2 and is covered by a path Z : [0, 1] — > M 4 connecting the 
points 1(0) = b\ and 1(1) = b . Thus either cot c*l (* denoting concatenation 
of paths) is a loop in M 4 based at b , and in both cases the image w.r.t. 
n M4 ^2 is freely homotopic to c. This shows that the image subgroup Im 
equals (n^i*), (tti(M 4 )) = vri(T 2 ) and thus (H). 

We need the above information in order to construct fast orbits in a given 
asymptotic direction. But we also need orbits with asymptotic speed zero in 
order to construct our 'stop and go orbits'. 

To that aim we notice that the kernel Ker((n A .f 4jT ) !|t ) C 7Ti(M 4 ) of our 
homomorphism is non-trivial. Indeed it contains the commutator subgroup 

Comm(7r 1 (M 4 )) = {ghg~ x h' x \ g, h G 7Ti(M 4 )} 

of the fundamental group, which is the smallest normal subgroup F with 
7Ti(M 4 )/F abelian. The fundamental group of M 4 is non-abelian, since the 
genus Q(M 4 ) > 3. So there exists a nontrivial 

s e Ker((n M4 ,T) J , s ^ id. 

By shortening a loop in the conjugacy class of s G 7Ti(M 4 ), we obtain a closed 
geodesic s : S 1 — > M 4 , our stop geodesic. After a reparametrization of time 
the projected geodesic IIm4,t ° s : S 1 — > T is a solution curve of the flow $* 
with initial conditions x s G and period T s . This stop orbit has asymptotic 
velocity v(x s ) = 0. 

Now in order to show our denseness result, we consider an arbitrary ve- 
locity i)6l 2 with modulus 

.. „ ^ V2(E - y max ) 
\\v\\ < — 

and seek for any e > an x G with ||u(xo) — ^|| < e - 

This is easy if v = because then we set xq := x s for our stop orbit with 
v(x s ) = 0. So we assume that v ^ and first approximate the direction 
v/\\v\\. 

26 



(16) 



We can find a lattice vector £' G C = tti(T) with 



IK'I 



e 

< - 



\[E — V me 
2 V2(E - V m 



(17) 



Given £', we now construct a closed geodesic g : S 1 — > M 4 whose associated 
periodic orbit on T starting at some point i 9 6 Eg has asymptotic direction 



v(x B 



\v(x 



and speed 



\v(x c 



> 



V2(E - V w 



(18) 



(19) 



\f~E ^mean 

As in the proof of Thm. ^]l|.0 we consider the closed straight lines 

k-.S^T, k(r) :=q + T-£(modC) 

with direction I := 4£' and initial point q G T. Since the loop k is at least 
four-periodic, we can lift k to M 4 obtaining a loop k : S 1 — ► M 4 (namely, the 
lift of to T 2 is at least two-periodic and the branched covering IT M4 >T 2 is 
only two-sheeted). 

Similar to (|7|), by an appropriate choice of the initial point q we can 
ensure that the length L(k) = \\£\\ ■ f \J E — V{k{r))dr of the corresponding 
loop in the Jacobi metric is bounded by 



L(k) < 



E — V m 



By shortening the loop k we obtain a geodesic g : S 1 — > M 4 which projects 
to a closed orbit on the torus starting at some x g G Tie- 

By the argument already used in the proof of Thm. |3.1| .[^ the period T g 
of that orbit is < ||£|| J^E-vr™) so ^ e as y m Pt°^ c velocity v(x g ) = £/T g 
meets (P^|). (|T^) is immediate from the construction since £'/\\£'\\ = £/\\£\\. 

Now this go orbit is too fast for our purposes. Therefore we find integers 
p,g6N with 



Tg_ ( ||r(.r (/) 



1 



< 5 



(20) 
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and consider for n G N the group elements 



By curve shortening we find a geodesic o n : S 1 — > M 4 in the conjugacy class 
of o n . We denote the period of the unit speed reparametrized d n by T n , and 
the period of the corresponding closed orbit on the torus by T n , and claim 
that 



T 



lim — 

n^oo n(p -Ts + q- Tg) 

This follows from 



1. 



(21) 



1. the Anosov property of the geodesic flow on the unit tangent bundle 
T1M4 of {M^g^E) proven in \Y1\ and 



2. the formula T n J J' §dr [J' 2iL ._ V[ , l{r)j) 



dr 



for the period. So the 



time reparametrization factor 1/(2(E — V(q))), seen as a function on 
T1M4, is Holder continuous. 

By 1) the geodesic flow line of o n approximates the go geodesic exponentially 
in n, then switches to the stop geodesic in n-uniformly bounded time, ap- 
proximates that geodesic exponentially in n, and finally switches back to the 
go geodesic in n-uniformly bounded time. 

Thus by 2) the ratio of times in (|2~ID goes to one as n — > 00. 

Let x n G T>e be a point on the torus orbit corresponding to the stop and 
go geodesic o n . Then v(x n ) = nq£/T n and v(x g ) = £/T g so that by ( [Zip 

pT s 



lim 

n— >oo 



r 9 + (p/q)T s 
The choice (^0|) of p/q implies that 
\v(x n )\\ - \\v\\\ < e/2 



\v(x 



<j) 



1 + 



? T 9 



lim 

n— >oo 



(22) 



for <5 > small. 

The geometric inequality 



vix n 



together with (|I 
large. 



< |||v(^n)|| 

r 



+ 



VlXr. 



V[X r , 



and (p2|) gives the result ||f — t> (x n )|| < £ for n 

□ 
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5 Semiclassics: Smooth Potentials 



We now compare the quantum system in the semiclassical limit with the 
classical one and thus mimick the definitions of section 3. 

The Schrodinger operators H h (k) on L 2 (T), fceT* have the eigenvalues 
E^(k). The semiclassical asymptotic velocities are defined by 

-h(i\ / ^ -1 V 'kE%(k) gradient exists 
n ^ ' ' \ otherwise. 

We equip the semiclassical phase space P h := N x T* with the semiclassical 
measure \ h := (2nh) d fii x /i 2 , where ji\ denotes counting measure on N and 
fj,2 Haar measure on the Brillouin zone T*. 

In order to compare classical and semiclassical quantities, we introduce 
the energy-velocity map 

A h :P h ^ R d+1 with A h {n, k) := (E%(k) , v%(k)) 

and the image measure v n := \ h {A K Y x . 

Our conjecture, which we shall prove in some special cases is: 
Conjecture 5.1 For all C-periodic potentials V G C°°(lR d , R) 
w* — lim v h = v 

(which means 

hm I f{x)du\x) = 
h\o J R d+i 

for continuous functions f e Cg(IR d+1 ,IR) of compact support). 

As can be already seen from d = 1-dimensional case, the supports of the 
semiclassical measures v n are in general much larger than the one of v. If 
the bands do not touch, then R x {0} belongs to supp^) (since then by 
symmetry v^(0) = 0), whereas the classical motion is ballistic above V mSuX . 

In this section we draw conclusions from Birkhoff's Ergodic Theorem 
which for some potentials imply the truth of our conjecture. 



/ f(x)du(x) 
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We first show that the range of semiclassical asymptotic velocities is in- 
cluded in the convex hull of the classical ones. No assumption on the inte- 
grability or ergodicity of the classical system is made. 

This involves a limit T — > oo, h — > which is controlled by the Birkhoff 
type proposition |5\2 . 

Let X be a compact metric space, consider a continuous flow 

$* : X -> X (teR) 
and a continuous map 
: X -> M d . 

Denote by M(X) the set of Borel probability measures on X, M(X, $) C 
M(X) the set of flow invariant ones and 

O t (x) :=W Oo &{x)dt. (23) 

By Birkhoff 's Theorem the good set 

G := \x G X I lim Or fx) exist and are equal } 

has measure //(G) = 1 for all // G M(X, $). We set O := limr^oo Or f G . 

The limit dist(0 T (x), 0(G)) -T is in general not uniform in x. How- 
ever this is true for the convex hull 

conv (0(G)). 

Denote for C C R d and for e > by C £ C R d the ^-neighbourhood of C, then 
it holds: 

Proposition 5.2 For all e > £/iere exists T £ > suc/i £/iet£ 

Or(X) C conv(0(G)) £ (|T| > T e ). 

Proof of |5.2| . By compactness of X and thus of 0(X) we could otherwise 
find an e > 0, a sequence of points x n G X and of times T n with T n — > ±oo 
such that 

2 := lim Tn (x n ) G M d 

n^oo 
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exists and z G" C e . W.l.o.g. we assume that T n — > +00. 

Consider the sequence of probability measures //„ G M(X) given by 

ti n (U) := i- |0 G [0,T n ] I G 17} | (C/ C XBorel). 

J- n 

We now use the following facts (see Thm. 6.10 in Walters' book, p2|): M(X) 
and M(X, $) are non-empty, convex, and compact in the weak-*-topology. 
The extreme points of M(X, $) coincide with the ergodic measures. By going 
to a subsequence, if necessary, 

fx ■= w * - lim fj, n G M(X) 

exists by compactness of M(X). \i G M(X, $) and as OT n (x n ) = J Odfi n 
the expectation f x Odjj, = z G" C £ . By Choquet decomposition of \i we would 
find an ergodic measure v G M(X, $) with J x Odv ^ C e . On the other hand 
by ergodicity of v there exists an x G X with 0(a;) = J x Odu, which is a 
contradiction. □ 

For I C K compact the phase space region 

Pj := {x G P I G /} 

is compact and ^'-invariant so that we are in the situation of Prop. |5.2| . 
The semiclassical analog of the thickened energy shell Pj is 

Pf:={{m,k)eP h \E h m {k)eI}. 

We equip them with the probability measures 

\t := - — - — on Pj 

KPi) 

and (for h small) 

A?:=^-^ on PI 
\\Pf) 

These induce the image probability measures \ii := Xjv" 1 and ^ := Aj(-u ft ) _1 
on the space M. d of asymptotic velocities. 
We shall now consider for e > intervals 

I E := [E — e,E + e]. 

and show our semiclassical results on the group velocities: 
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Theorem 5.3 Let V G C°°(M d ,R) be C-periodic, E G R,e > 0. 

1. Let C := coiw(v(Pj 2e )) . Then for all rj > 0, a.e. fceT* 3h~o V7i < ^ 

t# G if £y (A;) G I 6 ; 

2. letS := conv(supp(yUf 2 J) C M. d be the convex hull of the support of fii 2s , 
then the semiclassical measures concentrate inside S: For allrj > 

lim/ij (5„) = 1. 



Remark 5.4 In general S G C is much smaller than C. As an example for 
ergodic motion one has by Thm. \3. 1\ $ S = {0}, whereas by Thm. \3. C 
contains a disk of radius ^^-v^) ^ 

V V mean 



Proof. The proof of 1. is based on Theorem |5.2| , whereas for 2. we use the a.e 
convergence to the asymptotic velocity and a Shnirelman type argument. We 
shall freely use the semiclassical calculus as exposed in pOfl , || and references 
therein. 

First we state a lemma about Bloch decomposition of Anti-Wick quanti- 
zation. 

Lemma 5.5 Let f G C 6 °°(R 2d , R), /(p, g + £) = /(p, g) (£ G £), 

cAW„/, /" /-/^ „u / / „;,\ ^P^*? / , r2/Tn,cf\ 



f AW 1>:= f(p,q)<j> m (M) if-rh £ L (Mr)) 

where <t> p , q (x) ■= e~^ pq e^ px (p(x — q), <j)(x) := (7rh)~ d ^ 4 e~^ . 
It holds 

Uf A Wu -l = [ e f^ {k) M 



with 



f AW {k)i,:= If&qW^kXU^k)^)^ 
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Proof. By periodicity of / and unitarity of U we have 

Uf AW ^ = W f(p, q)U<f> p , q+l {U<j> p , q+i , Uj>)£%. 

eeC Jp (2vrft) 

Now 

U<j> Ptq+i (k,x) = e^+^e-'^U^^x) 

so 

uf Aw m = 

5 ttc jT * 1 1 ( 27r ") 

the claim follows now from Fourier inversion and the £*-periodicity of 

k' ^ (U<j) m (k'),Uij(k')}. ' □ 

A corollary of this lemma, the Egorov Theorem and the Weyl-Anti-Wick 
correspondence is: 

e iH*Wh f AW(Q e -iH*W* = (/ o ¥) AW {k) + T (h). 

Denote A h £ {k) = {j- Ef(k) e Q and X & C °°(R, R), supp X C x t/ a = 
1. For j e A^(fc) consider the eigenfunction of H h (k) and its Husimi 
distribution : P — > R, 

pj fc (p,g) := (27r^)- d |([/0 Pi(? (fc),^) L2(T) | 2 . 
We then have (using notation (|23|)) 

vf(k) = jMH)p)Tdp* k + T (/i). 



We apply Prop. |5.2| with X = Pi 2e . By time reversal symmetry G 
conv(-y(X)) so we find a T such that 

( X (H)p) T (X) C conv^X))^ 

and an h such that Ot{K) <r)/2. Thus 1. is proven. 
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Now we show that for a.e. k, for all r] > 

,. # {m e Ag(fc); |^(fc) - /xCgjgpg^ < ??} 1 

nm 7, i fc / 7 \ ; = 1 (24) 

This implies 2. as J x{H)vdp^ k C conv(supp/// 2 J by absolute continuity 

of x(H)dpj k w.r.t. A/ 2e , by time reversal symmetry and the bound 1 for the 
fraction. 

By Birkhoff 's Theorem 



T-+oo 

So the set 



lim / \vt(x) — v(x) \ d\i 2e = 0. 



B(T,rj) := {x G P he | |t) T (£) - > |} 



of phase space points x eventually giving rise to vt{x) S v /2 can be made 
small: 

\ 2£ (B(T,r)))< S - (T>T(S)). 
On the other hand, by convergence on the ^-independent Pj e : 
w* - lim - - pt b = ^h- 

For h < h(S) and T >T(5) we thus have 
^ £ A,* (B(T,r,))<S. 

By Tchebycheff 's inequality 

# [m e A r M | p h ^ k (B(T, V )) > VS} < V~5 #A?(fc) . 
For m in the complementary set it holds: 

f p X{H){ PT - v)dp^ k <R + 2|| X (^)^||ooV / 5 <\ (6< 6(rj))). 

To summarize: for a > there is a T and a set G\ C such that for 

h small enough \dj — J x{H)vdp"j k \ < rj for j E and > 1 — a. This 
finishes the proof of 2. □ 
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Corollary 5.6 If the classical motion is non-ballistic with probability one on 
an energy interval I: \ii = 5q, then Conjecture \5. 1\ holds true: 



w 



lim v 



v. 



For example, this is the case if the classical motion is ergodic. 



6 Semiclassics: Separable Potentials 

If the potential is separable, the distribution of semiclassical group velocities 
converges rapidly to the classical velocity distribution. We begin with the 
case of one dimension, and thus consider the operator H := —h 2 -^ + V{x) 
with potential V G C r (IR, E), assuming w.l.o.g. that V(x + 1) = V(x). The 
band function of the ra-th band for the quasimomentum k G [— 7r, it] is de- 
noted by E n {k) = E^(k). Of course ^§jr{k) = or the band functions touch 
at the band edges k = or ±7r. However, apart from small neighbourhoods 
of these values of the quasimomentum it holds: 

Proposition 6.1 Assume that the periodic potential V G C r (IR, M), r > 2. 

1. Then all bands in the energy interval [V max + e, oo) meet the following 
uniform estimate. If the quasimomentum \k\ G [h( r ~ 1 ^ 2 , it — Tfi"* 1 ^' 2 }, 
then 



sign (A;) • fT 1 ^(k) - (-l) n v cl (E n (k)) 
dk 



for some c = c(e) > 0, 



< ch 



-i 



Vd(E) = \ (jj ^-i— cftj (E > y max ), 

v c i(E) — (E < Knax) being the absolute value of the classical velocity. 
2. In the energy range [V^n, V ma , x — e] and for all k G [— n, n] 



35 



3. w* - lim v h = v. (25) 

Proof. 

1) We first consider the energy interval [V^ax + e, oo). For the k values under 
consideration, the Bloch eigenf unctions have no zeroes. So we are looking 
for zero-free solutions ip : R — > C of the differential equation 

Hep = Ecp. (26) 

The complex phase S := 5ln(</?) of such a solution solves the differential 
equation 

{S') 2 - ihS" -W = (27) 
with W := E — V . We solve this equation, using the Ansatz 



S(x) = S r (x) + tf +1 R(x, K) With S r (x) := 

With 



X ) 

n=0 



S (x) := / y/W$)dt 
Jo 

the recursion equation 

n 

S^x - £ S^-, = (n = l,...,r) (28) 



z=o 

has the continuous solution 



g»W = |/; ' S --" ) - E j|| ; " )S --"' ^ (»=l,...,r), 



since ^ = y/W > 0. 



In particular we have Si(x) = iln(y W(x)) + c. We set c := 0. Then on 
bounded intervals 



<p(x) = exp (i £ y/WQjdt/h) + O (H 1 ) . 



As a consequence of ( 28|) S n is real if n is even and imaginary if iV is odd. 
S' is 1 -periodic. 
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But unlike the real part of S r the imaginary part is always 1-periodic. 
This can be seen, e.g., by considering the formal power series Soc in h. By 
( |2TD the real part R := 9ft (iS^,) of the derivative is related to the imaginary 
part / := ^(S'J) by 

HR' = 2IR, 

so that / = ^h(hx(R))'. Thus the formal power series 9ft (Soo) is 1-periodic, 
which implies the same periodicity for its coefficients Szn+i- 

Of course this argument is even valid if (by finite differentiability of V) 
only finitely many coefficients S n are defined. 

Thus for E > V max 

<p r := exp ^SV^ = A r exp (^U^j : M. — > C 

is a function with periodic modulus A r > and phase U r = 9ft (S r ). We 
compare this with the solution ip of (^) with the initial values 

<p(0) := r (O), V'(0) := #(0). 
On bounded spatial intervals one has the estimate uniform in x and E 

<p(x) = (p r {x) + O (h r ), tp\x) = Cp' T (x) + o (h r ). 
The same is true for the matrix 



M(x) 



<p(x) ip(x) 
(p'(x) (p'{x) 



of the corresponding fundamental system (ip is linearly independent of ip). 
The monodromy matrix T := M(l) • M(0) _1 has determinant 1 and trace 

Tr(T) = 23fr(0 r (l)/£ r (O)) + O (h r ) = 2 cos(77 r (l)/ft) + O {ff). 

Now we consider those quasimomenta k for which 

dist(C/ r (l)/n,7r-Z) > h (r '- 1)/2 . 

For them |Tr(T)| < 2 if h < Hq so that the quasiperiodic function x 1— ► <£>(x) 
is bounded. Thus is a Bloch function with quasimomentum k, 

cos(» = §Tr(T). 
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Differentiating both sides w.r.t. k yields 

sin(ife) = sm(U r (l)/H) ■ ■ r ^ + O (H r ). (29) 

On the other hand 



sin(fc) = y 1 - (|Tr(T)) 2 = ± sin(C/ r (l)/ft) • (1 + O (h 1 )), 
so that ( ^9|) implies the relation 

n-^^r^V'^W (30) 



dk \ dE 

for the group velocity. By U r (l) = ^W{t)dt + (9 (/j 1 ) we get the estimate 
for quasimomenta € [^ r ~ 1 ^ 2 , 7r — ft^ -1 " 2 ] . 

2) For E < V ma , x — e we have fcK-E 1 ) = 0. Quantum mechanically it is well- 
known that the wave function and its derivatives are exponentially decreasing 
w.r.t. h well inside the potential well, say, for V(x) < e/2. This implies 
exponential decay of the group velocity. 
Namely 

h' 1 -^^) = j jn(x)dx = j n (x) 

for the current j n (x) = —ih^s(Lp n (k)(x)-^(p n (k)(x)) of the eigenfunction (p n (k) 
with eigenvalue E n {k) < V max — e, since the divergence of the current of 
eigenfunction vanishes. 

If we evaluate j at x inside the potential well, then we see that it is 
exponentially small. For more precise estimates valid in the multidimensional 
case we refer to Outassourt TH 



3) The convergence of the semiclassical measures v h to v follows from the 
following reasoning. 
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By a Weyl estimate we have weak-*-convergence in energy distribution: 
w* - limA^^r 1 = XH- 1 

with E h : P h — ► R, E h (m, k) := E^(k) being the energy function on the 
Fermi surface. The first two parts of the proposition exclude the energy 
interval (V max — e, Vm ax + e). However, as H \ 0, we can let £ = e(7i) \ 0, 
too. By the above Weyl estimate we do not loose anything of v in the 
semiclassical limit 

lim V(2vrn)- 1 \{k G T* | E*(k) - V max | < e(h)}\ = 0. 

meN 

Then ( p^) follows from 1) and 2). □ 

Corollary 6.2 Let V G C 2 (M a! ,]R) be a separable periodic potential. Then 
Conjecture \5. 1\ holds true: 

w* — lim v h = v. 

h\o 

Proof. By our assumption the potential is of the form 

d 

3=1 

with Vj G C 2 (IR, E) of some period lj > 0. Let Uj and Vj denote the 
(semi) classical measures for the one-dimensional potential Vj. Then 

u = (u 1 ,...,u d )L- 1 and v n — (i/f , . . . , vf}L~ x 

for the linear map L : (M. 2 ) d — > R d+1 , 

L(hi,vx, h 2 , v 2 ,..., h d , v d ) (ht + . . . + h d , v u ... , v d ). 

Although for d > 1 the linear map L is not injective and thus not proper, its 
restriction to 

x^ =1 ([^ >in ,oo) xR) (31) 

has this property, so that preimages of compactly supported functions are 
still compactly supported. We may restrict L to (|3lD, since the support 
of (ui, . . . , v d ) is contained in (|3ij), the spectrum of Hj being contained in 

[Vj, m in, OO) 

So multidimensional convergence follows from the one-dimensional one. 

□ 
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